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Abstract 

The propagation of a pulse in a nonlinear array of oscillators is influenced 
by the nature of the array and by its coupling to a thermal environment. For 
example, in some arrays a pulse can be speeded up -while in others a pulse can 
be slowed do-wn by raising the temperature. We begin by sho-wing that an 
energy pulse (ID) or energy front (2D) travels more rapidly and remains more 
localized over greater distances in an isolated array (microcanonical) of hard 
springs than in a harmonic array ot in a soft-springed array. Increasing the 
pulse amplitude causes it to speed up in a hard chain, leaves the pulse speed 
unchanged in a harmonic system, and slo-ws do-wn the pulse in a soft chain. 
Connection of each site to a thermal environment (canonical) affects these 
results very differently in each type of array. In a hard chain the dissipative 
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forces slow down the pulse while raising the temperature speeds it up. In a 

soft chain the opposite occurs: the dissipative forces actually speed up the 

pulse while raising the temperature slows it down. In a harmonic chain neither 

dissipation nor temperature changes affect the pulse speed. These and other 

results are explained on the basis of the frequency vs energy relations in the 

various arrays. 

PACS numbers: 05.40.Ca, 05.45.Xt, 02.50.Ey, 63.20.Pw 
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I. INTRODUCTION 

In recent years there has been a great deal of interest in the interplay of nonlinearity 
and applied forcing (deterministic and/or stochastic) in the stationary and transport prop- 
erties of discrete spatially extended systems [|l|. The ability of discrete anharmonic arrays 
to localize and propagate energy in a persistent fashion, and the fact that noise may act 
(sometimes against one's intuition) to enhance these properties, has led to particularly in- 
tense activity [0-^. Interesting noise- induced phenomena include stochastic resonance [Q, 
noise-induced phase transitions |^, noise- induced front propagation [0, and array-enhanced 
stochastic resonance p. 

Our interests in this area have been motivated by the relative dearth of information 
concerning the effects of a thermal environment on the sometimes exquisite balances that 
are required to achieve these interesting resonances and persistences |]3|, p^JTT[] . At the same 



time, we have also noted that most of the literature has concentrated on overdamped ar- 
rays (often motivated by mathematical or computational constraints rather than physical 
considerations), a restriction that leaves out important inertial effects and that is easily 
overcome. 

Perhaps the simplest generic discrete arrays in which to analyze these issues are sys- 
tems of oscillators consisting of masses that may be subject to local monostable potentials 
(harmonic or anharmonic) and nearest neighbor monostable interactions (harmonic or an- 
harmonic) (other generic arrays of current interest are bistable units linearly or nonlinearly 
connected to one another). These are the systems of choice in our work, and we have sep- 
arated our inquiries into three distinct groups of questions: 1) The study of such arrays 
in thermal equilibrium [jlO|. The questions here concern the spatial and temporal "energy 
landscape" that determines the degree of spontaneous energy localization due to thermal 
fluctuations and the temporal persistence of high or low energy regions; 2) The study of the 
propagation of a persistent signal applied at one end of the array [0]. The questions here 
concern the signal-to-noise ratio and distance of signal propagation; 3) The study of the 



propagation of an initial (5- function energy pulse (this work). The questions here concern 
the velocity of propagation and the dispersion of such a pulse. 

It is useful and relevant to provide a very brief summary of our conclusions on the first 
two sets of questions. Our work on equilibrium energy landscapes ^^ was based on chains 
of harmonically coupled oscillators subject to a local potential that may be anharmonic. 
Each oscillator is connected to a heat bath at temperature T. We analyzed the thermal 
fluctuations and their persistence as influenced by the local potential (we compared hard, 
harmonic, and soft potentials), the strength of the harmonic coupling between the oscilla- 
tors, the strength of the dissipative force connecting each mass to the heat bath, and the 
temperature. Among our conclusions are the following: 1) An increase in temperature in 
weakly coupled soft chains leads not only to greater energy fluctuations but also to a slower 
decay of these fluctuations; 2) An increase in temperature in weakly dissipative hard chains 
leads not only to greater energy fluctuations but also to a slower decay of these fluctuations; 
3) High-energy-fluctuation mobility in harmonically coupled nonlinear chains in thermal 
equilibrium does not occur beyond that which is observed in a completely harmonic chain. 

However, we noted earlier that interest in energy localization in perfect arrays, as con- 
trasted with localization induced by disorder, arises in part because localized energy in these 
systems may be mobile. Dispersionless or very slowly dispersive mobility would make it pos- 
sible for localized energy to reach a predetermined location where it can participate in a 
physical or chemical event. Our results raised the possibility of observing such localized mo- 
bility if the anharmonicity lies in the interoscillator interactions rather than (or in addition 
to) the local potentials. We ascertained that a persistent sinusoidal force applied to one 
site of a chain of masses connected by anharmonic springs may indeed propagate along the 



chain |Tl[]. Furthermore, we demonstrated a set of resonance phenomena that we have called 
thermal resonances because they involve optimization via temperature control. In particu- 
lar, these results establish the existence of optimal finite temperatures for the enhancement 
of the signal-to-noise ratio at any site along the chain, and of an optimal temperature for 
maximal distance of propagation along the chain. These resonances differ from the usual 



noise-enhanced propagation where the noise is external and/or the system is overdamped. 

This work addresses the third set of questions posed above concerning the way in which 
a nonequihbrium initial condition in the form of an energy pulse propagates as the system 
relaxes toward equilibrium. More specifically, we investigate the motion and dispersion of 
such an energy pulse and the effects of finite temperatures on pulse propagation. In view of 
our earlier results on thermal resonances, perhaps the most interesting question to be asked 
at this point is this: Is it possible to enhance pulse propagation via temperature control? 

In order to monitor the evolution of the nonequihbrium initial condition it is useful to 
partition the Hamiltonian as 

H = j:En (1) 

n 

where En contains the kinetic energy of site n and an appropriate portion of the potential en- 
ergy of interaction with its nearest neighbors (1/2 in one dimension, 1/4 in two dimensions). 
In one dimension 

n^ 1 1 

En = Y + ^V{Xn+l,Xn) + -V{Xn, X„_i) . (2) 

In Section ^ the potentials considered in this paper are briefly presented. Section ^ 
contains our analysis and results for one-dimensional oscillator chains. Here we discuss ways 
to characterize the mobility and dispersion of an initial localized impulse, and compare 



the behaviors of harmonic, hard anharmonic, and soft anharmonic chains. In Section |I^ we 
present some results for isolated two-dimensional arrays and note some interesting geometric 
features with perhaps unanticipated consequences. Section is a summary of results. 

II. POTENTIALS 

The particular potentials as a function of the relative displacement y = Xn — Xn-i used 
in our presentations are the harmonic, 

Vo{y) = \y\ (3) 



a hard anharmonic, 

VH{y) = \y\ (4) 

and a soft anharmonic, 

Vs{y) = k[\y\-Hl + \y\)] . (5) 

The anharmonic potentials have been chosen to be strictly hardening and strictly softening, 
respectively, with increasing amplitude. The potentials are shown in the first panel in Fig. |l|. 
In almost all our simulations we take k = 1. 

The displacement variable ^ of a single oscillator of energy E m a, potential V{y) satisfies 
the equation of motion 

dy 



^^ ±p[E-V{y)]. (6) 

This equation can be integrated and, in particular, one can express the period of oscillation 
t{E) and the frequency of oscillation uj{E) as 

^{E) Jo ^2[E-V{y)] 

The amplitude ymax is the positive solution of the equation V{y) = E. The resulting 
oscillation frequencies obtained from the integration of Eq. (|^) for the three potentials with 
k = 1 as well as that of the frequently used "quadratic plus quartic potential" are shown in 
the second panel of Fig. |l| fl^ . 
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FIG. 1. First panel: the potentials defined in Eqs. (^)-@ with k = 1. Sohd curve : harmonic 
potential, Vb(y). Dashed curve: hard anharmonic potential, Vh{y)- Dotted curve: soft anharmonic 
potential, Vs{y). Second panel: frequency as a function of the oscillator energy for these potentials. 
The dot-dashed line shows the frequency curve for the commonly used potential Vo(y) + Vh{y)- 



The frequency vs energy variations seen in Fig. |l| can be shown via rescaling and bounding 
arguments to represent general features of hardening and softening monostable potentials. 
The exercise is trivial if the potential is of the form 



since then 



whence 



t{E) = 4 



Viy) = -/ 
a 



dy 



(8) 



-1 



l/a .1 



dz 



E-^ 



B^E^-^ 



a/2 1 



(9) 



/ ^N Svr ^1 1 
uj{E) = -^E2-^ 



(10) 



The coefficient Ba can be expressed exactly in terms of the beta function and is equal to 2% 
for the harmonic potential. 



If the potential is not of the simple single-power form it is still possible to bound the 
resulting energy dependence to establish the trend [0. For example, the soft potential @ 
is bounded below by {k/2)\y\ and above by k\y\. These bounds immediately lead to the 
conclusion that the associated uj{E) must decrease as E^^^"^. The argument for a mixed 
power potential such as V{y) = ^y^ + ^y^ is a bit more cumbersome but otherwise similar: 
by making the change of variables from y to ||/^ + ^y* = Ez'^ one can show not only that 
uj{E) is an increasing function of E but that it lies above the harmonic potential result for 
any positive E. 

Figure |l] summarizes the well known frequency characteristics of oscillators: for a har- 
monic oscillator the frequency is independent of energy (and, with our parameters, equal 
to unity); for a hard oscillator the frequency increases with energy, while that of a soft 
oscillator decreases with energy. The hard oscillator frequency curve starts below the other 
two if a harmonic portion is not included. These frequency-energy trends are generalized 
to oscillator chains in Appendix ^. The frequency vs energy behavior will figure promi- 
nently in our subsequent interpretations. In particular, the following broad view seems to 
be overarchingly supported: the speed and dispersion of pulse propagation in discrete arrays 
of oscillators are principally dependent on the mean frequency associated with the energy 
in the pulse. Higher frequencies lead to faster propagation and slower dispersion. 

III. ONE-DIMENSIONAL ARRAYS 

We consider one-dimensional arrays of 2N+1 sites numbered from — A^ to N with periodic 
boundary conditions. We distinguish isolated chains (that is, ones not connected to a heat 
bath), chains connected to a heat bath at zero temperature, and finite temperature chains. 
This provides an opportunity to organize the effects of different parameters on the behavior 
of the chains. 

In all cases at time t = a kinetic energy e is imparted to one particular oscillator (the 
oscillator at n = 0) of the chain. If the chain is isolated or at zero temperature, this initial 
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impulse is applied to an otherwise quiescent chain. At finite temperatures the chain is first 
allowed to equilibrate and then this impulse is imparted in addition to the thermal motions 
already present. We then observe how this initial 5- function impulse propagates and spreads 
along the chain, and how these behaviors depend on system parameters. 

A. Isolated Chains 

The equations of motion for an isolated chain are 

d 

X'n = -^ [V{Xn - X„_i) + 1^(X„+1 - Xn)] ■ (H) 

The initial conditions are 

2^n(0) = for all n, 

i;„(0) = forn ^ 0, 

xo{0) =po = V2^ . (12) 

For a harmonic array this system can of course be solved exactly, and we do so in Appendix 0. 
The analytic harmonic results are helpful and informative, although our discussion is primar- 
ily based on simulation results since the anharmonic chains cannot be solved analytically. 
The numerical integration of the equations of motion is performed using the second order 
Heun's method (which is equivalent to a second order Runge Kutta integration) ||T3| , p^ with 
a time step At = 0.0001. 

One can think of the dynamics ensuing from the initial momentum impulse in two equiv- 
alent ways. One is to interpret the a;„ and i;„ as displacements and momenta along the 
chain. Two symmetric pulses start from site zero and move to the left and to the right along 
the chain, and our discussion focuses on either of these two identical pulses. This symmetry 
occurs regardless of the sign of the initial momentum since the energy does not depend on 
the sign, i. e., the contraction of the spring between sites n = and n = 1 that follows 
an initial positive impulse has exactly the same effect as the equal extension of the spring 



between sites n = and n = —1. Alternatively, one can think of x„ and Xn as displacements 
and momenta perpendicular to the chain, the sign then simply representing motion "up" or 
"down" . The symmetry around the site n = is then even more obvious. 

In any case, the energy e excites the displacements as well as momenta of other oscillators 
as it moves and disperses. The evolution can be characterized in a number of ways. We have 
found the most useful to be the mean distance of the pulse from the initial site, defined as 



(X)^ 



2-m ^n 



(13) 



and the dispersion 



2 _ / 2\ / \2 I^n'^ En , ,2 

a = (x ) — {x) = —— — [X) 



(14) 



(The sums over n extend from —N to N .) Here the E^ are the local energies defined in 
Eq. (^ and, since these depend on time, so do the mean distance and the variance. The 
time dependence of the mean distance traveled is a measure of the velocity of the pulse, and 
that of the dispersion is a measure of how long the pulse survives before it degrades to a 
uniform distribution. An indication of the progression of a pulse is shown in Appendix ^ 
for a harmonic chain. 
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FIG. 2. First panel: Mean distance traveled by the initial energy pulse as a function of time 
for the hard, harmonic, and soft chains with several initial momentum amplitudes. Second panel: 
pulse dispersion as a function of mean pulse position for the three chains with initial amplitude 
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Results for the mean distance traveled by the pulse as a function of time for isolated 
chains of 151 sites are shown in the first panel in Fig. ^ for the hard, harmonic, and soft 
potentials and for various values of the initial pulse amplitude po- The mean distance varies 
essentially linearly with time in all cases (this is only approximately true in all cases - 
even the harmonic oscillator exhibits early deviations from linear behavior due to inertial 
effects). The important results apparent from Fig. |^ are summarized as follows and can be 
understood from the frequency vs energy trends in Fig. ||: 

1. The pulse velocity in the harmonic chain is independent of the initial amplitude. This 
reflects the energy-independence of the mean frequency (and in fact of the entire 
frequency spectrum) for harmonic chains (also see Appendix B). 

2. The pulse velocity in the hard chain increases with increasing initial amplitude. This 
is because the mean frequency for the hard chains increases with increasing energy. 

3. The pulse velocity in the soft chain decreases with increasing initial amplitude. This 
is because the mean frequency for the soft chains decreases with increasing energy. 

We note that with our choice of potentials the velocity in the hard chain for very weak 
initial amplitudes may actually lie below that of the harmonic chain or even the soft chain 
because we have omitted a harmonic contribution to the hard potential, but the hard chain 
velocity necessarily increases and surpasses that of the other chains with increasing initial 
pulse amplitude. 

Not only is the pulse transmitted more rapidly in the hard isolated chains than in the 
others, but the pulse retains its integrity over longer distances in the hard chain. This is 
seen in the second panel in Fig. ^ The dispersion a^ is shown for the three chains for a 
particular initial pulse amplitude. Rather than the dispersion as a function of time, the 
dispersion is shown as a function of position along the chain so that the pulse widths at a 
particular location along the chain can be compared directly. Clearly the hard chain pulse is 
the most compact at a given distance from the initially disturbed site (a plot of a^ vs t would 
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show the opposite trend, that is, the pulse in the hard chain would have the greatest width, 
but it will have traveled a much greater distance than the pulses in the other chains). This 
combination of results leads to interesting geometrical consequences in higher dimensions 
(see Sec. 0). 

B. Chains at Zero Temperature 

If the chains are connected to a heat bath at zero temperature, the equations of motion 
Eq. (|ll]) are modified by the inclusion of the dissipative contribution, 

d 



Xyi 



dXn 



[V{xn - a;„„i) + l^(x„+i - Xr, 



/•^n 



(15) 



where 7 is the dissipation parameter. The initial conditions are as set forth in Eqs. 
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FIG. 3. Mean pulse position as a function of time for the hard, harmonic, and soft chains with 
initial amplitude pq = 8, with and without friction. 

The mean distance traveled by the pulse is shown in Fig. |^ for each of the chains with 
and without friction so that the frictional effects can be clearly established. The salient 
results can again be understood from the frequency vs energy trends in Fig. |1|: 
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1. The pulse velocity in the harmonic chain is independent of friction. This again reflects 
the energy-independence of the mean frequency for harmonic chains. The energy loss 
suffered through the frictional effects therefore does not affect the pulse velocity. 

2. The pulse in the hard chain slows down with time in the presence of a frictional force. 
This is because the chain loses energy via friction, and the mean chain frequency 
decreases with decreasing energy. 

3. The pulse in the soft chain speeds up with time in the presence of a frictional force. This 
is because the chain loses energy via friction, and the mean chain frequency increases 
with decreasing energy. 

The dependence of the pulse width on friction (not shown explicitly) follows trends that 
are consistent with our other results. An increase in friction causes the pulse to narrow in 
the soft chain. This is consistent with the observation that higher frequencies are associated 
with narrower pulses. In a harmonic chain there is also some narrowing of the pulse, but not 
nearly as much as in the soft chain (detailed explanation of this would require consideration 
of the spectrum beyond just the mean frequency). In the hard chain we can not make 
an unequivocal claim from our numerical results because the dependence of pulse width 
on friction for our parameters is extremely weak, with perhaps a very small amount of 
narrowing. 

C. Chains at Finite Temperature 

If the chains are connected to a heat bath at temperature T, the equations of motion 
Eq. (|15D are further modified by the inclusion of the fiuctuating contribution, 

d 

X'n = --K [^(a^n - Xn-l) + V^(x„+i - Xn)] " 'JXn + r/n(t) . (16) 



^n 



The rjnit) are mutually uncorrelated zero-centered Gaussian (5-correlated fluctuations that 
satisfy the fluctuation-dissipation relation: 
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(77„(t)) = 0, {Vn{t)vAt')) = 2^kBTSnjS{t - 1') . (17) 



The initial conditions are now no longer given by Eqs. (|T2|). Instead, the chain is allowed 
to equilibrate at temperature T and at time t = an additional impulse of amplitude \/2e 
is added to the thermal velocity of site n = 0. The integration of the equations of motion 
proceeds as before, but now we report averages over 100 realizations. The system is initially 
allowed to relax over enough iterations to insure thermal equilibrium, after which we take 
our "measurements." 

The pulse dynamics is no longer conveniently characterized by the mean pulse velocity 
(although this was the most useful and direct characterization in the absence of thermal 
fluctuations). This is because there is now a thermal background that causes fluctuations 
and distortions of the information in this mean (as well as in other simple moments and 
measures such as the pulse maximum). We flnd that the most suggestive presentation of 
the dynamics is that of the energy proflle itself. An illustrative set of typical proflles for 
chains of 51 sites is presented in Fig. ^, showing energy proflles as a function of time on 
the S*'^ site on either side of n = as a function of temperature. In all cases there is a 
delay time until the pulse reaches the S*'^ site (reflecting a flnite velocity). The local energy 
around this site then reaches a maximum, and the pulse moves on, leaving behind a series 
of later energy oscillations at ever decreasing amplitudes that eventually settle down to the 
appropriate thermal levels. The after-oscillations are derived analytically in Appendix ^ for 
the harmonic case. The discussion below concentrates exclusively on the flrst pulse, which 
we think of as characterizing the arrival of the disturbance at that site. 
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FIG. 4. Energy profile vs time at the 5* site for hard and harmonic chains at different tem- 
peratures. The damping parameter in all cases is 7 = 0.2 and the initial pulse amplitude is po = 8. 

The important conclusions, some illustrated in the figure, can once again be understood 
from the trends in Fig. |1| and include the following: 

1. The pulse velocity in the harmonic chain is independent of temperature. This is illus- 
trated in the figure by the fact that the peak of the pulse reaches the particular site 
under observation at the same time for the two temperatures shown. The reason once 
again is that the characteristic frequencies of the chain are independent of energy and 
therefore the inclusion of thermal effects is immaterial to this measure. 

2. The pulse velocity in the hard chain increases with increasing temperature. This is 
illustrated by the ever earlier arrival of the pulse at the site under observation with 
increasing temperature. The reason is that the mean frequency of the chain increases 
with energy, so that the hard chain at higher temperatures is associated with a higher 
frequency than at lower temperatures and hence with a faster pulse. 

3. The pulse velocity in the soft chain decreases with increasing temperature. This is not 
explicitly illustrated in the figure, but is due to the decrease of the mean frequency 
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with increasing temperature. Thus the soft chain at higher temperatures is associated 
with a lower frequency and hence a slower pulse. 

4. The hard chain not only transmits pulses more rapidly than the other chains, increas- 
ingly so with increasing temperature, but it also transmits the most compact and 
persistent pulses at any temperature. This is seen not only by the obviously smaller 
width of the pulses in the hard chain, but by the fact that the energy trace "left be- 
hind" as the first pulse passes through is lower in the hard chain than in the other 
cases. 

The pulses in all cases become more dispersive with increasing temperature. This be- 
havior is clearly evident in Fig. ^ for the hard and harmonic chains, as is the fact that the 
temperature dependence of the pulse width is weakest for the hard chain (and strongest 
for the soft chain). These dependences complement those described earlier for the pulse 
width as a function of friction: increasing friction in all cases narrows the pulse (subject to 
our caveat concerning the hard chain mentioned earlier) while increasing the temperature 
broadens it, both of these dependences being weakest for the hard chain. 

IV. TWO-DIMENSIONAL ISOLATED ARRAYS 



We showed in Sec. |IIIA| that a pulse travels more rapidly and less dispersively in an 
isolated hard chain than in a harmonic or soft chain. In higher dimensions these two ten- 
dencies, that of moving faster and that of maintaining the energy localized, leads to some 
interesting geometric effects and to very different pulse propagation properties depending 
on the spatial configuration of the initial condition. 

In one dimension one could visualize the displacements and momenta x, x as describing 
motion along the chain or perpendicular to the chain. In two dimensions these are distinct 
cases: a generalization of the first requires introduction of two-dimensional coordinates 
{x,y) and momenta {x,y). The second requires only a single perpendicular coordinate z 
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and associated momentum i for each site, and this is the case we pursue. We thus consider 
a two-dimensional square array of dimension (2A^ + 1) x (2A^ + 1) wherein motion occurs in 
a direction perpendicular to the array. The Hamiltonian with Znj = Pn,j is expressed as a 
sum of local energy contributions, 

H = Y.En,, (18) 

where 



.2 

_ t 



771 _ Pn,j ^ 



k' 



\.^n,j ^n+l,j ) ~r \^n,j ^n—l,j) ~r \^n,j ^n,j+l) T \^n,j ^n,j—l) 



\Zfij Zfi^ijj -f- yZfij Zfi—ijj -f- \Zfij ZnJ+l) 1 \^n,j ^n,j — l) 



(19) 

For the harmonic case we take k = 0.5 and k' = 0, and for the hard anharmonic array 
we set A; = and k' = 0.5. Our lattices are of size 51 x 51 and our integration time 
step is At = 0.0005. The boundary conditions are immaterial (although we happen to 
use boundaries whose edge sites have only three connections and its corners sites only two) 
because the lattices are sufficiently large for the initial excitations not to reach the boundaries 
within the time of our computations. 

We consider two initial excitation geometries. In one a "front" is created by exciting all 
sites along the line (0,j), — A^ < J < A^, with the same initial momentum poj = Po- The 
front then moves symmetrically away from this line and its motion is measured by the mean 
distance and dispersion (in all double sums n and j each range from — A^ to A^) 

(-) - %^^ , (20) 

2 / 2\ / \2 2^n,j^ ^n,] 



a^ = {x^) - {xY = ^^^-r^ - {xY . (21) 

In the other, an initial pulse of kinetic energy is deposited at the central site of the array. 
We then measure the mean radial distance of the pulse from the origin, 
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r = 



En,,- V^^ + fEnJ 



(22) 



(the dispersion in this case is less informative but can also be monitored if desired). The 
motion and dispersion in this geometry are expected to be roughly spherically symmetric 
subject to the square connectivity of the lattice. 
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FIG. 5. Snapshot at a subsequent time of the energy distribution for the propagation of an 
initial front. Upper panel: harmonic lattice. Lower panel: hard lattice. 





FIG. 6. Snapshot at a subsequent time of the energy distribution for the propagation of an 
initial pulse at the center of the array. First panel: harmonic lattice. Second panel: hard lattice. 

Typical gray-scale snapshots of the energy distribution are shown in Figs. |^ and |^, and 
the differences, while easily understood, are clearly dramatic. In the case of the front, the 
tendency of a hard lattice to propagate faster than the harmonic lattice while maintaining 
the energy more localized is clearly realized. The associated mean distance and dispersion 



that quantify the comparison are shown in Fig. ^ In the case of an initial point pulse, on 
the other hand, there is clearly a conflict between rapid motion and smaller dispersion - one 
can be realized only at the expense of the other. The latter "wins" : the pulse remains more 
localized in time in the hard lattice than in the harmonic. The associated mean radius is 
shown in Fig. ^. In the anharmonic lattice the pulse at first expands as fast as the harmonic 
but it essentially quickly saturates while the harmonic pulse continues to disperse. 
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FIG. 7. First panel: temporal evolution of the mean distance (x) of front propagation in 
harmonic (solid line) and hard anharmonic (dashed line) lattices. Second panel: the associated 
dispersion as a function of the mean distance traveled. 




FIG. 8. Temporal evolution of the mean radius (r) of pulse propagation in harmonic (solid line) 
and hard anharmonic (dashed line) lattices. 
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V. CONCLUSIONS 

In this paper we have considered pulse propagation in discrete arrays of masses connected 
by harmonic or anharmonic springs. We have focused on the pulse velocity and width, and 
have found a pattern of behavior that can be strongly correlated with the energy dependence 
of the mean array frequency. 

First we investigated the propagation of pulses in isolated (microcanonical) arrays. We 
found that in a hard array an amplitude increase causes a pulse to travel more rapidly and 
less dispersively. In a harmonic array the pulse speed and width are independent of pulse 
amplitude, while in a soft array a more intense pulse travels more slowly and spreads out 
more rapidly. These trends are a result of the fact that in a hard array the mean frequency 
increases with energy, in a harmonic array it is independent of energy, and in a soft array 
the mean frequency decreases with increasing energy. In higher dimensions these trends lead 
to interesting initial condition dependences that in turn may lead to apparently "opposite" 
behavior in different cases. Thus, for example, a front in a two-dimensional isolated hard 
array propagates more rapidly and more sharply than in harmonic or soft arrays, and the 
effect is enhanced if the front is more intense. On the other hand, a point pulse in a hard 
array spreads more slowly than in the others: it is not possible in this geometry to both 
propagate quickly and yet retain a strong localization of energy, and the latter tendency 
dominates the dynamics. 

We then investigated the effects on pulse propagation of connecting the nonlinear chains 
to a heat bath (we did this only for the ID arrays). We found that dissipative forces tend to 
slow down the pulse in the hard array, leave its speed unchanged in the harmonic chain, and 
actually speed it up in the soft array. This somewhat counterintuitive behavior is, however, 
fully consistent with the observation that dissipation causes a decrease in energy and hence 
a decrease in mean frequency in the hard case and an increase in mean frequency in the soft 
chain (and no change in the mean frequency of the harmonic chain). Dissipation in all cases 
causes a narrowing of the pulse, the effect being greatest in the soft array. 
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An increase in temperature has the opposite (and again at first sight perhaps somewhat 
counterintuitive) effect: it speeds up the pulse in the hard array, leaves it unchanged in the 
harmonic array, and slows it down in the soft chain. Again this behavior is consistent with 
the frequency vs energy trends and the fact that an increase in temperature is associated 
with an increase in the energy of the chain. A temperature increase in all cases causes a 
broadening of the pulse, the effect again being greatest in the soft array. 
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APPENDIX A: FREQUENCY VS ENERGY FOR OSCILLATOR CHAINS 

Consider a chain of oscillators, and let us focus on the displacement variable x of a 
particular mass in the chain, say oscillator j, whose displacement satisfies the equation of 
motion 

1/2 

dt 



E -^V{Xn-Xn^l] 






The period of oscillation for oscillator j can be defined in analogy with Eq. (p^: 
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where x' stands for the set of all the x's except Xj, and similarly for p'. The upper limit of 
integration Xmax depends not only on E but on all the other displacements and momenta, 
and is the positive value of Xj at which the denominator of the integrand vanishes. The 
resulting u with all the coordinate and momentum dependences is not very useful, but it 
would seem reasonable to simply average over all possible values of these coordinates and 
momenta and thus obtain an average period. We define the average period as 

r(i?) = (r(E;x',p')) 

/ ■ ■ ■ / rfx7 ■ ■ ■ / dp- /„— dXj (2 [E - En V(Xn - Xn_i)] - En^j Pn) "'^^ 

/■■■/rfx'/---/dp' 

(A3) 

The limits of integration not explicitly indicated are appropriate nested relations among the 
variables and the energy such that the argument of the square root always remains positive. 
The multiple integral in the denominator covers the same integration regime and insures 
proper normalization for this average. Our interest lies in extracting the energy dependence 
- the remaining energy-independent coefficients are complicated and not important for our 
arguments. If the pair potentials are powers as in the single oscillator example, the scaling 
argument can be generalized by introducing scaled variables Zn ~ {xn — a;„_i)i?^/" and 
Un ~ PnE^^'^ with appropriate constants of proportionality. The limits of integration then 
become independent of energy and the only energy dependence arises from factoring an 
^1/2 from the square root in the denominator and an E^/°' from the numerator because it 
contains one z— integration more than the denominator. The result, as before, is that 

t{E) = B^E^--^ (A4) 

with a complicated but energy-independent expression for the coefficient Ba, and therefore 

^(E) ^ ^- E^^^" ■ (A5) 
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More complicated potentials require suitable generalization of this argument, but the 
result in any case is that the average frequencies for the hard, harmonic, and soft chains 
follow the same trends as those shown in Fig. |1|. 

APPENDIX B: ISOLATED LINEAR OSCILLATOR CHAIN 



Although linear oscillator chains are of course fully understood, it is nevertheless useful 
to present aspects of their behavior in the context of the present discussion. 

The linear equations of motion (|TT]) with the initial conditions ( [T^ ) are easily solved: 
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where the frequencies Uq obey the dispersion relation 



uj„ = 4k sin 



2 f 27rg 
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V2A^+L 

and where Jn{z) denotes the Bessel function of the first kind of integer order n. (The energy- 
independence of the frequencies for the harmonic chain seen here in the e-independence of 
the Uq is prominent in our discussions throughout this paper.) The momenta are then 

(B3) 



Xn{t) = V^ J2ni2Vkt). 
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FIG. 9. Temporal evolution of the energy En{t) at several sites in a harmonic chain with force 
constant k = 1 and initial momentum po = 2. 

Using a number of relations obeyed by the Bessel functions it is possible to combine these 
results and obtain for the local energy the simple expression 



Enit) 



jl{2Vkt) + \ji+,{2^t) + \ji-ii:^Vkt) 



(B4) 



The energy profiles for various sites are shown in Fig. ^. The n = 5 profile (here obtained 
from the analytic expression (|B4| ) also appears in Fig. ^ (there obtained by numerical in- 
tegration). Note that the energy is not transported in a single absorption-emission process 
but rather in a series of oscillatory steps of decreasing amplitude. Our analysis in the body 
of the paper focuses on the first energy pulse. 
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FIG. 10. First panel: Mean distance traveled by the pulse (dashed curve) and pulse maximum 
(solid curve) as a function of time for a harmonic chain of unit force constant. Second panel: Pulse 
maximum as a function of time for k = 2 (long-dashed curve), k = 1 (solid curve), and k = 1/2 
(short-dashed curve). The slopes of these numerically generated curves are essentially as given in 
the analytic expression Eq. (p6|). 



In Section |111A| we rely on (x(t)), the mean distance traveled by the pulse as a function of 
time, as one measure to characterize the transport properties of our arrays. An alternative 
measure that can be calculated analytically for the harmonic chain (but turns out to be 
somewhat less convenient for numerical computation) is the time-dependent site n*{t) at 
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which the energy is a maximum. Because the passing energy pulse in general leaves a track 
behind it, one expects n*{t) = Xmax{t) to grow more rapidly than {x{t)). That this is indeed 
the case is illustrated in the first panel of Fig. [10|, where both quantities are shown for a 
harmonic chain with unit force constant. The steps in the Xmax curve are a consequence 
of the discreteness of the problem. The analytic result for n*{t) is obtained by maximizing 
Eq. ( p4D with respect to n and is, after some manipulation, found to be the solution of the 
relation 



l + VJTT ^ J2ni2Vkt) 

4:71 2v^tJ2„_i(2v^t) ■ 

Except for a very short initial transient the solution is essentially linear in time and exceed- 
ingly simple: 

n*{t) = Xraaxit) ^ Vkt . (B6) 



This dependence is confirmed in the second panel of Fig. ^ for three values of the force 
constant. The curves shown are obtained numerically, and differ from the analytic straight 
lines only at the very earliest times by an exceedingly small barely visible amount. 
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